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The  reactants'  and  products’  wavefunctions  are  calculated  from. their 
respective  distortion  potentials  as  one-term  adiabatically  separable 
approximations.  Both  the  distortion  potentials  and  the  interaction  poten- 
Jtials  are  extracted  straightforwardly  from  the  LEPS  surface.  The  novel 
feature  of  the  approach  is  that  for  the  first  time  accurate  results  for  the 
absolute  values  of  the  reaction  probability  are  obtained  from  a  simple 
overlap  of  single-channel  approximate  wavefunctions  obtained  directly 
from  the  respective  parts  of  the  potential  energy  surface  for  the  reac- 


1.  INTRODUCTION 


The  rates  of  many  chemical  reactions  at  low  and  moderate  tempera¬ 
tures  are  determined  mainly  by  the  reactive  transition  probabilities  of 
-energies  below  and  around  the  classical  reaction  threshold.  At  such  ener¬ 
gies  the  scattering  is  predominantly  inelastic  and  the  reactive  scattering 
probabilities  are  low.  One  can  expect  therefore  that  treating  the  reactive 
scattering  as  a  perturbation  to  the  much  simpler  inelastic  scattering 
may  yield  a  satisfactory  simple  method  for  evaluating  the  reactive  transi¬ 
tion  probabilities  in  the  threshold  region,  and  hence  for  obtaining  ther¬ 
mal  reaction  rates. 

Treatments  of  this  type,  known  as  Distorted  Wave  Born  Approxima¬ 
tions  (DWBA)  for  reactive  collisions,  have  been  extensively  studied  in 
recent  years.1*5  They  consist  in  evaluating  the  reactive  transition  proba¬ 
bility  as  an  integral  over  two  separate  nonreactive  wavefunctions,  one 
corresponding  to  the  reactants'  and  one  to  the  products’  configuration, 
and  an  exchange  interaction  operator.  Further  simplified  versions  of  the 
treatment,  known  as  Franck-Condon  theories  of  chemical  reactions,  have 
also  been  frequently  considered®*17  One  simplification  consists  in  treating 
the  interaction  operator  as  a  constant,  which  reduces  the  expression  to 
an  overlap  integral  between  the  two  distorted  wavefunctions.  Utilizing 
analytical  approximations  to  the  distorted  wavefunctions,  usually  in 
terms  of  Airy  functions  for  the  unbound  motion  and  harmonic  oscilator 
wavefunctions  for  the  bound  motion,  also  allows  calculating  an  approxi¬ 
mation  to  the  transition  probability  in  closed  form,  clearly  at  a  further 
loss  of  accuracy. 

At  energies  below  the  classical  reaction  threshold  the  main  contribu¬ 
tion  to  the  DWBA  integral  comes  from  the  exponentially  decaying  tail  end 


of  the  distorted  wavefunctions!8  As  a  result,  small  displacements  of  the 
tail  position  can  lead  to  errors  of  several  orders  of  magnitude  in  evaluat¬ 
ing  the  DWBA  integral?’11  Most  of  the  work  done  using  the  DWBA  approx¬ 
imation  in  reactive  scattering  has  been  actually  concerned  with  evaluating 
the  vibrational  distributions  of  the  products,  and  the  absolute  values  of 
the  transition  probabilities  have  only  rarely8,5,12  been  evaluated. 

More  recently!9'23  it  has  been  shown  that  the  DWBA  formalism  can 
yield  accurate  absolute  values  for  the  reactive  transition  probabilities  in 
the  threshold  region  if  properly  evaluated  distorted  wavefunctions  and 
interaction  potential  are  employed  for  the  reactants  and  the  products. 
We  have  developed  a  simple  DWBA  treatment  of  degenerate  and 
nearly  degenerate  H-atom  transfer  reactions  in  which  the  distorted 
wavefunctions  are  constructed  as  a  linear  combination  of  two  adiabati- 
cally  separable  wavefunction.  The  treatment  has  been  shown  to  give 
excellent  results  for  several  model  systems!9'21  Hubbard,  Shi  and  Miller23 
have  demonstrated  that  the  DWBA  approximation  can  give  accurate 
results  in  the  threshold  region  for  the  H  +  //2  reaction,  if  converged  dis¬ 
torted  wavefunctions,  extracted  from  an  accurate  coupled  channel 
numerical  solution  of  the  complete  Schrodinger  equation  for  the  prob¬ 
lem,  and  the  exact  interaction  operator,  are  utilized  in  the  evaluation  of 
the  integral. 

In  this  paper  we  develop  a  simpler  version  of  our  earlier  DWBA  treat¬ 
ment  of  the  light  atom  transfer  reactions!9'21  The  present  work  is  closer 
in  spirit  and  simplicity  to  the  Franck-Condon  treatments6-17  but  still 
retains  the  high  accuracy  of  the  transition  probabilities  for  the  light 
atom  transfer  reactions  mentioned  above!9*22  It  is  also  not  restricted  to 
transfer  between  two  degenerate  or  nearly  degenerate  vibrational  states. 
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The  distorted  wavefunctions  are  constructed  as  single  channel  adiabati- 
cally  separable  eigenfunctions  of  readily  defined  distortion  potentials. 

The  basic  features  of  the  present  model  are  given  in  Sec.  II.  Numeri- 
-cal  tests  of  the  treatment  for  model  potential  energy  surfaces  for  the  I  4 
HI  and  Cl  4  HBr  reactions  are  given  in  Sec.  III.  The  results  are  discussed 
in  Sec.  IV. 


H.  DIABATIC  FORMALISM 

We  consider  a  rearrangement  collision  which  involves  a  transition 
from  a  vibrational  state  n  of  the  reactant  molecule  BC  to  a  vibrational 
state  m  of  the  product  molecule  AB, 

A  4  BC(n)  -*  AB(m)  4  C  .  (2.1) 

The  Schrodinger  equation  can  be  conveniently  written  in  terms  of  the  two 
alternative  sets  of  Jacobian  mass-weighted  coordinates  [rp,Rp\ 

appropriate  for  describing  the  reactants'  and  the  products’ 
configurations: 

H  =  Ha  -  -  i--^r-4  ha(ra,Ra)  ;  A„(r „,/?„)  =  -  ^--~4  Va{ra,Ra) 


2  dRl 


2  0r‘ 


H  =  -  j-  -g-4  hp{rp,Rfi)  ;  hfi{rfi,Rp)  =  -  -■  ^-4  Vp{rp,Rp)  . 


2  dRj 


where  rt  =  ri',  rt'  being  the  internal  coordinate  of  the  diatomic 
molecule  (i  =  a  for  reactants,  0  for  products);  m*  is  the  reduced  mass  of 
the  molecule  i.  Similarly  R^  =  M^RX\  Rx'  being  the  separation  between 
the  atom  and  the  center  of  mass  of  the  diatomic  molecule  i;  Ux  is  the 
reduced  mass  of  the  relevant  atom-diatom  pair,  hifo,/?*)  is  a  vibrational 
Hamiltonian  for  fixed  Ri . 


v.vv.  . 


V.W  ' V  • vs 


In  the  asymptotic  region  the  dependence  of  the  vibrational  Hamil¬ 
tonian  hi(rxRx)  on  /?,  vanishes.  The  two  h,'s  reduce  to  the  vibrational 
Hamiltonian  for  a  single  diatomic  molecule: 

hBc(ra)  =  -  Y  +  v-(r«)  ■  h*B (rfi)  =  ~  Vfrp)  (2.3) 


2  6r 


In  the  case  of  transfer  of  a  light  particle  lying  between  two  heavy 
ones  (i.e.,  mB  «  mA,mc)  the  polar  coordinates  (p,0 ), 

P  =  U?f  +rf)»  (2.4) 

0  =  Lan~l(ra/  Ra)  =  em  -  tan-,(r^//?^)  ,  (2.5) 

are,  as  shown  elsewher^,9-21,24  particularly  well-suited  for  simplifying  the 
approximate  treatment  of  the  dynamics,  due  to  the  near-separability  of 
the  motion  along  the  fast  (0)  and  slow  (p)  coordinates.  The  angle  0m 
between  the  skewed  axes  in  Fig.  1  equals  tan“,[m^/(maAfa)}*].  The  Ham¬ 


iltonian  in  these  polar  coordinates  is  written  as 

H  _  _  i__af _ 1 _ d _ 1  92 

2  dp2  2P  dp  2 p2  d0z 


•+  V(0.p) 


It  is  also  convenient  to  define  an  arc  length,  measured  from  the  Ra  axis 
in  Fig.  1,  and  to  denote  the  maximum  arc  length  by  sm: 


s-p0,  sm=p0m  (2.7) 

Asymptotically  the  system  can  be  in  the  arrangement  channel  a  or  /B. 
For  large  values  of  p  one  can  relate  the  asymptotic  form  of  Eq.  (2.6)  to 
that  of  Eq.  (2.2)  by  employing  the  approximate  relations  between  the 
polar  and  the  Cartesian  coordinates  in  the  region  relevant  to  the  motion 
of  the  system,  namely  in  the  regions  s  «  sm  and  sm  -  s  «  sm.  We  have 
there 
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Ra « P  .  ra  w  s  ;  Rp*tp  ,  rfi*sm-s  .  .  (2.0) 


In  order  to  obtain  the  DWBA  reaction  amplitude  for  the  rearrange¬ 
ment  reaction  we  divide  the  potential  energy  function 

V(s.p)  =  Va(s,p)  =  V„(s,p)  ,  (2.9) 

into  two  parts 


1  ?.p)=  V'(s,p)  +  V"(s,p) 

V^.p)=Vl(s,p)  +  Vj/(s,p) 


(210) 


The  distortion  potential  V{(s  ,p)  includes  all  the  distortion  of  the  potential 
between  the  atoms  in  the  diatomic  molecule  due  to  the  proximity  of  the 
third  particle  but  none  of  the  interactions  responsible  for  the  rearrange¬ 
ment.  The  interaction  potential  V(!{s  ,p)  on  the  other  hand  represents 
solely  the  reactive  part  of  the  interaction  which  causes  the  rearrange¬ 
ment  (t  =  a,/?). 

We  take  the  nonreactive  wavefunctions  for  both  arrangements  to  be 
the  solutions  of  the  Schro  dinger  equations  with  the  distortion  potential 
v£  only. 
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(211) 


In  the  absence  of  the  reactive  interactions  the  motion  along  p  and  along  s 
is  nearly  adiabatically  separable  within  each  of  the  valleys  (reactants’, 
products’)  of  the  potential  energy  surface.  Thus,  a  good  approximation  to 
these  nonreactive  solutions  of  Eq.  (2.11)  can  be  obtained  using  the  adia¬ 
batic  approximation,4 


-  8- 


1 ^n(s-p)s  <finip)Xn(s  ,p) 

^(s,p)a  <p£i(p)x&(s  ,p)  . 


(212) 


where  the  x’s  are  the  eigenfunctions  of  the  nonreactive  Hamiltonian  (Eq. 
(2. 1 1))  for  a  fixed  p  and  depend  parametrically  on  p.  They  satisfy 


-  |rri"+  Xn(s-p)  =  Cn(p)x“(s.p) 

2  ds* 

-  vfi(s>P)  X&(s.p)  =  e&(p)xA(s.p)  • 

2  as*1 


(2.13) 


where  we  have  used  ds  =  pdfi  in  the  derivative.  The  distorted  single 
channel  vibrational  wavefunctions  x(s  ,p)  are  functions  of  s  which  depend 
parametrically  on  p.  It  can  be  readily  shown  that  the  p-dependent 
coefficients  tp*  in  Eq  (2.12)  satisfy  approximately  the  Born- 
Oppenheimer-hke  Eq.  (2.14). 


1  a2 

2  dp2 

1  a2 

2  dp2 


Vm(p)  =  0 


(2.14) 


The  DWBA  transition  amplitude  for  the  reaction  is  then  given  by 

Tf:l  =  <1r&\VH\i,Z>  ,  (2.15) 

where  V'n(V'm)  is  the  wavefunction  for  the  reactants  (products)  when  the 
potential  is  the  nonreactive  part  v£(v£).  A  perturbative  treatment  that 
includes  the  reactive  part  of  the  potential  V11  to  the  first  order  can  be 
expected  to  give  accurate  results  at  energies  below  the  classical  thres¬ 
hold,  where  the  reaction  can  be  viewed  as  a  perturbation  to  the  nonreac¬ 
tive  scattering. 


Using  the  single  channel  adiabatic  approximation  (Eq.  (2. 13)).the  T- 
matrix  expression  (Eq.  (2.15))  for  going  from  state  n  of  channel  a  to  state 
m  of  channel  fi  can  be  rewritten  as 

■  »•  *m 

T$1  =  f  dp  <P&(p)<P%(p)f  ds  Xm(s  >p)  Va  (s  ’P)Xn(s  >p)  •  (2.16) 

0  0 

For  every  value  of  the  radial  coordinate  p  between  the  asymptotic 
region  and  the  saddle  point  region,  the  cut  along  the  s  coordinate  on  the 
potential  energy  surface  of  a  light  atom  transfer  reaction  represents  a 
double  minimum  potential.  We  confine  our  attention  to  such  p’s  for  which 
the  cut  is  a  double  minimum  potential.  The  boundary  between  "reac¬ 
tants"  and  "products"  at  any  finite  p  is  defined  as  the  top  of  the  barrier 
s  =  s6(p)  in  Fig.  2.  We  now  define  V £  in  such  a  way  as  to  be  identical  to 
Va  in  the  reactants’  valley  and  to  be  a  constant,  namely  the  value  of  Va  at 
s6  (p)  in  the  products’  valley.  Similarly,  Vlp  is  defined  at  each  p  to  be  Vp 
in  the  products'  valley  and  to  be  Vp[sb  ,p)  in  the  reactants’  valley.  That  is, 


vi  (s  )  _  v«(s  -P)  s  <  sj  (p)  t 
a(  Va(s*,p)  s  >s*(p)  • 


(2.17) 


Similarly  for  Vp(s,p )  we  have 


Vm(s  o)  =  Vs -P)  s^s*(p) 
*  ,P  Vp(sb,p)  s  <  sb(p) 


(2.18) 


The  reactive  part  of  the  potential  V{l(p,s)  is  then  given  by 

K//=  Vt  -  V[  (i=a,/J) 


(2.19) 


V»  vanishes  for  s  s=  sb(p),  for  all  p,  and  Vjj  vanishes  for  s  sb(p)  for  all 
p.  Equation  (2.16)  then  reduces  to 


TiSfi  =  /  dp  <ph(p)<p%(p)  Jf  ds  xfh(s  ,p)V'J (s  ,p)xn(s  .p)  .  (2.20) 

0  »r( ft) 
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UI.  APPLICATION  TO  I  +  HI  AND  Cl  +  HBr  SYSTEMS 

The  major  part  of  the  contribution  to  the  integral  above  comes  from 
a  narrow  range  of  values  of  p  near  the  classical  turning  point.  We  shall 
-therefore  approximate  the  solution  of  Eq.  (2.14)  by  the  solution  of  the 
equivalent  equation  in  which  the  potential  ( E  -  ej),  where  i  =  a,/5  and 
j  =  n,m,  is  linearized  around  the  classical  turning-point  (p°  and  p£  for 
the  reactants’  and  products’  channels,  respectively): 

(E  -  e“)  =  (p  -  pZ)Fn 

(3.1) 

(E  —  e£)  =  (p  —  Pm)Fm 

The  Fj ’s  are  the  derivatives  of  E  —  Cj  with  respect  to  p  at  the  turning 
point  p?.  That  is,19 

=  (2/  Bnj*  Ai[-Bn(p  -  p°)] 


(32) 


(2 / Bm)*  Ai[-Bm(p  -  p^)]  . 
where  Bj  =  (2 Fj)z/3  and  Ai  is  the  Airy  function. 
The  reaction  probability  is  given  by 
p£  =  [2tt 


T£&  =  p-8^2— w- /  dp  Ai [ -Bm (p  ~ Pm)]Ai[-Bn{p-p^)} 

x /  —x5.(s.p)[Va(s,p)  -  Vg  (s,p)]x"(s.p)  ■ 

«*  p 

The  integral  over  s  in  Eq  (3.4), 

V$i(p)  =  f  ,p)[Va(s  ,p)  -  Vla  (s  ,p)]Xn(s  >P)  ds/p 


(3.3) 


(34) 


(3.5) 


was  evaluated  numerically.  The  eigenfunctions  Xn  and  xln  of  Eq.  .(2.13) 
were  also  evaluated  numerically  using  a  finite  element  method.  Since 
VX&  is  approximately  an  exponential  function  of  p  in  the  region  around 
JLhe  turning  points,  a  region  making  the  principal  contribution  to  the 
integral,  it  can  be  written  as 

KSS  3  •  (3  6) 

where  is  the  value  of  V£ £(p)  at  the  "average”  turning  point  p#  half¬ 
way  between  the  turning  points  p§  and  p£.  The  slopes  Fm  and  Fn  for  the 
present  system  are  almost  identical 

Fm~Fn  =  F;  Bm*Bn=B  (3.7) 

Eq.  (3.4)  can  then  be  written  as 

T^  =  ^-fAi[B(p-pS)]e^-’,S)Ai[B{p-pg)]dp  .  (3.8) 

which  can  readily  be  evaluated  analytically19  to  yield 

=  p  exp  ,24 F  exp  "&F  •  (3'9) 

Here,  A  is  the  value  of  V“(p)  -  V*(p)  at  p§.  More  details  on  the  (3.8) 

integral  evaluation  are  given  in  Ref.  19,  Sec.  V.  In  the  case  of  symmetric 

H-atom  exchange  between  the  same  vibrational  states  the  same  formula 

applies  with  A  =  0  and,  then  Eq.  (3.7)  is  exact.  In  the  case  when  states  n 

and  m  are  far  from  being  degenerate,  Eq.  (3.4)  can  also  be  readily 
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evaluated  analytically  using  a  slightly  different  procedure 

The  forward  and  the  reverse  transition  probabilities  and 

obtained  from  Eq.  (3.0)  are  slightly  different  due  to  the  lack  of  symmetry 
in  the  approximate  Eq.  (3.3).  The  mean  of  the  forward  and  reverse  proba- 


bilitics  and  was  found  to  give  a  more  accurate  approximation  to 
the  transition  probability  than  either  P^m  or  Pmk  alone  and  is  used  below 
in  the  comparison  of  the  results.  Due  to  the  closeness  of  the  P£&.  and 
PjmL  results  the  particular  form  of  the  mean  used  (geometric  or  arith¬ 
metic)  did  not  influence  the  result  significantly. 

Equation  (3.8)  was  used  to  calculate  the  probabilities  for  H-atom 
transfer  in  two  different  collinear  systems:  I  +  HI  -*  IH  +  I  and 
CL  +  HBr  -*  CLH  +  Br  for  which  accurate  numerical  solutions  of  the  two- 
state  problem  are  available. 

The  transition  probabilities  obtained  by  using  Eq.  (3.8)  are  presented 
in  Figs.  3  and  4  as  a  function  of  the  total  energy  for  the  I  HI  and  ClHffr 
systems,  respectively.  The  results  of  an  accurate  numerical  coupled 
channel  solution  of  the  CLHBr  problem  are  also  given  for  comparison  in 
Fig.  4.  The  results  for  the  IHI  system  in  Fig.  5  are  compared  with  the 
results  of  an  accurate  numerical  solution  of  the  two  state  problem.19  The 
latter  have  been  shown  '  to  be  in  excellent  agreement  with  the  coupled 
channel  results  on  a  slightly  different  potential  energy  surface  for  the 
reaction  for  which  accurate  coupled  channel  calculations  have  been  per- 
formed  The  numerical  results  obtained  for  Figs.  3  and  4  using  (3.6)  were 
indistingishable,  in  the  Figures,  from  those  obtained  numerically  using 
(2.16). 

IV.  DISCUSSION 

The  results  presented  in  the  preceding  section  demonstrate  that  the 
present  simplified  DWBA  treatment  leads  to  very  accurate  results  for  the 
reactive  transition  probabilities  below  the  classical  threshold.  As  dis¬ 
cussed  in  Ref.  20  such  accuracy  in  the  threshold  region  is  sufficient  to 
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detcrmine  the  reaction  rates  at  low  and  normal  temperatures  for  the 
reactions  tested  as  well  as  to  explain  curvature  of  the  Arrhenius  plots  for 
H-transfer  reaction  at  low  temperatures15,16,20 

The  present  treatment  is  conceptually  simpler  than  the  one  we  dev¬ 
ised  in  Refs.  19-21  since  it  follows  a  very  straightforward  formulation  of 
the  DWBA  distorted  wavefunctions  and  interaction  potentials.  It  is  also 
computationally  simpler  since  it  does  not  involve  evaluation  of  eigenfunc¬ 
tions  of  double  minimum  potentials.  Although  the  eigenfunctions  of  the 
one-dimensional  anharmonic  potentials  are  evaluated  in  the  present  work 
numerically,  approximate  analytical  procedures  for  their  evaluation  can 
also  be  devised. 

The  validity  of  the  present  method  is  restricted  to  light  atom 
transfer  reactions  with  moderate  to  high  potential  energy  barriers  to  the 
reaction.  The  present  formalism  hinges  on  the  possibility  of  defining 
separate  vibrational  states  for  the  reactants  and  the  products,  and  can 
not  be  implemented  if  the  barrier  is  not  appreciably  larger  than  the 
zero-point  vibrational  energy  of  the  reactants.  The  systems  to  which  the 
present  formalism  does  apply  are  the  ones  for  which  the  tunneling  contri¬ 
bution  to  the  reaction  rate  is  significant  and  for  which  the  DWBA  transi¬ 
tion  probabilities  around  the  reaction  threshold  determine  the  thermal 
reaction  rates. 

The  formalism  presented  here  is  not  restricted  to  transfer  between 
nearly  degenerate  vibrational  states,  since  no  such  assumption  is  used  in 
the  derivation.  Comparisons  with  accurate  numerical  results  are,  how¬ 
ever,  presented  only  for  nearly  degenerate  and  degenerate  systems  due 
to  the  paucity  or  absence  of  accurate  results  for  other  systems. 
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F1GURE  CAPTIONS 
Fig.  1. 

Plot  indicating  the  various  variables  used  on  the  text.  The  ’skewing 
angle*  6m  is  also  given. 

Fig.  2. 

Profile  of  the  potential  energy  surface  for  a  fixed  p  is  depicted  by  the 
curve  ABCD.  The  nonreactive  surfaces  V, (s  p)  and  v£(s,p)  are  given 
by  the  curves  ABCG  and  FCDE,  respectively. 

Fig.  3. 

Log  pl!o  versus  total  energy  for  the  /  +  ///-»///  +  /  system.  The  cir¬ 
cles  denote  the  results  given  by  Eq.  (3.3)  and  the  solid  line 
represents  the  results  of  the  numerical  solution  of  the  two-state 
problem. 

Fig.  4. 

Log  Pqz  versus  total  energy  for  the  CL  +  HBr  -*  CLH  +  Br  system. 
The  circles  denote  the  results  given  by  Eq.  (3.3)  and  the  solid  line 
represents  the  results  of  an  accurate  coupled  channel  numerical 
solution  of  the  problem. 
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